Evidence for crossed Andreev reflection in superconductor-ferromagnet hybrid 

structures 
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We have measured the non-local resistance of aluminum-iron spin-valve structures fabricated by 
e-beam lithography and shadow evaporation. The sample geometry consists of an aluminum bar 
with two or more ferromagnetic wires forming point contacts to the aluminum at varying distances 
from each other. In the normal state of aluminum, we observe a spin-valve signal which allows 
us to control the relative orientation of the magnetizations of the ferromagnetic contacts. In the 
superconducting state, at low temperatures and excitation voltages well below the gap, we observe 
a spin-dependent non-local resistance which decays on a smaller length scale than the normal-state 
spin-valve signal. The sign, magnitude and decay length of this signal is consistent with predictions 
made for crossed Andreev reflection (CAR). 
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Singlet superconductivity and ferromagnetism are 
competing long-range orders which favor a different align- 
ment of electron spins, antiparallel and parallel, respec- 
tively. Therefore, they generally exclude each other in 
homogenous bulk materials. In mesoscopic hybrid struc- 
tures, the interplay of superconductivity and ferromag- 
netism leads to rich novel physics. Recent experimen- 
tal studies on equilibrium properties of superconductor- 
ferromagnet (SF) interfaces have shown that the local 
density of states (LDOS) on the superconducting side 
is strongly affected by the pair breaking effect of the 
ferromagnetic exchange field Q, and on the ferromag- 
netic side an oscillatory behavior of the LDOS due to 
the exchange splitting of the spin sub-bands has been 
observed [2j. Transport properties of SF point contacts 
0- El and FSF planar junctions 0] show a suppres- 
sion of spin-polarized current injection into the super- 
conductor, and in SFS Josephson junctions 0,0, Eli a 7r- 
state with a spontaneous equilibrium Josephson current 
is observed. In this work, we report on spin-dependent 
transport properties of the superconducting condensate 
on length scales comparable to the coherence length. We 
have combined ferromagnet-superconductor point con- 
tacts with a non-local spin-valve geometry, and find ev- 
idence for crossed Andreev reflection (CAR), i.e. the 
splitting of a Cooper pair into two spatially separated 
leads. 

Figure ^ shows an SEM image of one of our sam- 
ples, together with a schematic view of the experiment. 
The samples were fabricated by e-beam lithography and 
shadow evaporation techniques. First, 20 nm of iron is 
evaporated onto an oxidized silicon substrate to form a 
series of wires (vertical), with varying width (50-120 nm) 
and a small tip at the end. Their elongated shape creates 
a magnetic shape anisotropy which confines their magne- 




Fe 



u 



± 0- 




FIG. 1: (Color online) SEM image of sample T2, and ex- 
perimental scheme. Three vertical iron wires are connected 
by point contacts to a horizontal aluminum bar. The outer 
two iron wires have additional voltage probes. An example of 
the current injection (7a) and both local (J7a) and non-local 
(Ub) voltage detection scheme for one pair of contacts is also 
shown. 



tizations to be aligned along the wires, pointing either up- 
wards or downwards, with different coercive fields due to 
the width variation. In a second evaporation step under a 
different angle, and without breaking the vacuum, an alu- 
minum bar of 80 nm height and 250 nm width (horizon- 
tal), is created. It slightly touches the iron wires, forming 
metallic contacts of about 20 x 50 nm, much smaller than 
the dirty-limit coherence length £ ss 200 — 300 nm of the 
aluminum. This is crucial to avoid a suppression of su- 
perconductivity by the proximity effect. A current I a is 
injected through one of the contacts into the supercon- 
ductor, and a voltage Ub is detected by a second contact 
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outside the current path, defining the non-local resistance 
R = Ub/Ia- Samples of two different layouts were inves- 
tigated. Layout T (shown) has three iron wires and addi- 
tional voltage probes for detection of the local resistance 
Ua/Ia- Layout S has six iron wires over a larger range 
of distances, and no additional voltage probes. For each 
sample, different combinations of contacts have been used 
as injector-detector pair in order to study the dependence 
of R on the distance d between the two contacts. Here, 
we show data from one sample of each layout, T2 and S5. 
The samples were mounted into a shielded box thermally 
anchored to the mixing chamber of a dilution refrigera- 
tor. The measurement lines were fed through a series 
of filters to eliminate RF and microwave radiation from 
the shielded box. Resistance was measured with a low- 
frequency AC resistance bridge, with an RMS excitation 
amplitude 7 exc . For differential resistance measurements, 
an additional DC bias current was applied. 

For each pair of contacts, we first characterize the spin- 
valve behavior in the normal state. The inset of Fig. 
shows the non-local resistance for two contacts as a func- 
tion of the magnetic field B applied along the direction of 
the ferromagnetic wires at a temperature T — 1.8 K, with 
an excitation amplitude I cxc — 1 /iA. For each direction 
of the magnetic field, the data show two sharp stepwise 
changes of the resistance which correspond to the reversal 
of the magnetizations of the two contacts at their indi- 
vidual coercive fields, yielding each two states for parallel 
and antiparallel magnetization alignment. These states 
remain stable at zero external magnetic field. The resis- 
tance difference Ai?N between parallel and antiparallel 
alignment in the normal state is governed by the injection 
of non-equilibrium magnetization into aluminum, and the 
associated splitting of the chemical potentials for spin up 
and down electrons |l0| . A quantitative analysis of Ai?N 
will be given below. 

For each of the four magnetization states of a given 
contact pair, the sample is cooled through the super- 
conducting transition at zero external magnetic field, 
and subsequently warmed up again. Here a very small 
excitation amplitude / cxc = 50 nA is chosen to avoid 
both quasiparticle injection above the superconducting 
gap and self-heating. At this small excitation (and cor- 
respondingly small detector voltage amplitude Ub ~ 
0.5 nV), the limited absolute accuracy of the resistance 
bridge yields an artificial offset of w 150 mil to R. The 
observed data depend only on the relative magnetization 
alignment of the two Fe wires, and are therefore aver- 
aged over parallel and antiparallel configurations. The 
results for one pair of contacts is shown in Fig. |21 Above 
the superconducting transition temperature, the signal 
is constant, with the well resolved difference Ai?N be- 
tween parallel and antiparallel alignment. At the super- 
conducting transition at T c « 1.15 K, a large peak is 
observed, which is essentially independent of the magne- 
tization alignment. This peak is due to charge imbalance 
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FIG. 2: (Color online) Non-local resistance R as a function 
of the temperature T at zero external magnetic field for two 
contacts of sample T2 at a distance d = 310 nm. The solid and 
dashed line corresponds to parallel and antiparallel alignment 
of the injector and detector magnetizations, respectively. 
Inset: Non-local resistance R at T = 1.8 K as a function 
of a magnetic field B applied parallel to the iron wires, for 
two contacts of sample S5 at a distance d — 210 nm. The 
resistance jumps correspond to magnetization reversals of the 
ferromagnetic contacts. The arrows indicate the four different 
magnetization states. 

created by the injection of quasiparticles above the su- 
perconducting gap 0, 0| , which is still small close to 
T c . As the gap opens further towards lower tempera- 
tures, the charge-imbalance peak subsides, and finally, 
the resistance becomes constant again at T < 250 mK. 
In this regime, the signal for parallel alignment is still 
somewhat larger than for the antiparallel case. The dif- 
ference ARg between the non-local resistances for par- 
allel and antiparallel alignment in the superconducting 
state is the central result of this work, and we will now 
further clarify its physical origin. 

A possible source of the spin-dependent signal in 
the superconducting state are quasiparticles propagating 
above the superconducting gap. Such a signal has been 
observed in the charge-imbalance region near T c using 
current injection into a niobium film through a tunnel 
junction |l3j. It was found to be an order of magnitude 
smaller than the normal-state spin- valve signal at a tem- 
perature 0.3% below T c , and became unobservably small 
about 1.5% below T c . We do not observe such a signal 
close to T c , presumably because it is buried in the steep 
spin-independent increase caused by charge imbalance. 

To check further whether AR$ observed at low tem- 
perature is due to quasiparticles above the gap, or due 
to sub-gap transport, we have employed the additional 
voltage probes of sample T2 to characterize the point 
contacts. The local resistance Ua/Ia of one contact is 
shown in Fig. [31 as a function of temperature. The data 
can be described by R{T) = R + Rbtk(T) + R ci (T), 
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FIG. 3: (Color online) Symbols: Local resistance Ua/Ia vs. 
temperature T for a contact of sample T2. Solid line: fit to 
Ro + Rbtk(T) + R C1 (T) (see text). Dash-dotted line: the con- 
tribution R + Rbtk(T). 

Inset: Differential local resistance (UIa/cLIa vs. 7a for the 
same contact (symbols), and differential non-local resistance 
dUB/dlA for a detector at distance d = 210 nm (solid line: 
parallel, dashed line: antiparallel alignment). The dash- 
dotted line is the fit Ro + Rbtk{Ua) (see text), the dotted 
line indicates Ro + R pc in both frames. 

where Rq = 63.7 is the resistance of the ferromag- 
netic strip leading to the point-contact, Rbtk(T) is the 
temperature dependence of the point-contact resistance 
according to a spin-polarized extension of the BTK the- 
ory 0, and Rn@) is a contact resistance due to 
charge imbalance [lj, [la, LL21 UM > which is essentially the 
same as observed in the non-local resistance. Parameters 
for the BTK model are R pc = 6.2 O, P = 0.45, Z = 0.22 
and A = 177 //eV for the normal state point-contact re- 
sistance, the interfacial spin polarization, the interface 
transparency, and the gap, respectively. The fit to our 
model is excellent except for a slight broadening around 
T c . 

The low bias part (|/a| < 1-7 ftA) of the differen- 
tial local resistance dU^/dlx can be well described by 
Ro + Rbtk(Ua) using U A = Ia{Ro + R P c) (i.e. assum- 
ing insignificant inelastic scattering in the Fe strip) and 
the same parameters as for the temperature dependence, 
whereas the differential non-local resistance dU^/dl^ re- 
mains constant in that bias region (inset of Fig. |3J). At 
higher bias |/a| > 1.7 fiA, both local and non-local dif- 
ferential resistance show a steep increase due to charge 
imbalance, i.e. quasiparticle injection above the gap. The 
apparent reduction of the gap compared to the BTK pre- 
diction can be attributed to self-heating at higher bias. 
We conclude that the constant Ai?s observed in the low- 
bias region is due to sub-gap transport. 

Sub-gap transport between normal metals and super- 
conductors is mediated by Andreev reflection [l^: An 
electron from the normal metal enters the superconduc- 



FIG. 4: (Color online) Difference A_Rn,s between parallel and 
antiparallel alignment vs. distance d between the contacts, for 
two samples in both the normal (open symbols) and super- 
conducting (closed symbols) state. The solid and dashed lines 
are fits to Q and JjjJ, respectively, as described in the text. 

tor, and a hole of opposite spin is retro-reflected, cre- 
ating a spin-singlet Cooper pair in the superconductor. 
In a multi-terminal structure, it has been predicted that 
the incident electron and the retro-reflected hole may be 
transmitted through different contacts (crossed Andreev 
reflection, or CAR), as long as the distance between the 
two contacts does not exceed the superconducting co- 
herence length [2fl . EH | , yielding a negative non-local re- 
sistance. There is a second non-local process, electron 
cotunneling (EC), where an electron enters the super- 
conductor through one contact, and an electron of the 
same spin leaves through the second contact, yielding 
positive non-local resistance. If the two contacts are spin 
polarized, EC and CAR are favorable for parallel and 
antiparallel magnetization alignment, respectively. We 
propose the superposition of CAR and EC processes as 
the origin of the observed non-local resistance difference 
in the superconducting state. 

To obtain a quantitative description of our results, and 
to test our hypothesis, we have extracted Ai?N and Ai?g 
from the raw data and investigated their dependence on 
the distance d between the injector and detector contact. 
These data are shown in Fig. Both signals are of simi- 
lar magnitude at short distance, and decay at increasing 
d. The decay in the superconducting state takes place 
on a much shorter length scale than in the normal state, 
indicating the different transport mechanism. 

For the normal state data, we use a semiclassical model 
of spin injection and diffusion The key parameters 
are the bulk spin polarization p = (ct f ^— (t f ^)/((t f ^+(t f ^) 
of the ferromagnet, where tr F ^ and a F ^ are the spin- 
dependent conductivities of the ferromagnet, and the 
quantities K N = p N A^/A N and K F = p F X F /A\ where 
p( N > F ) aim As are the (spin-averaged) resistivity and 
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spin-diffusion length in the normal metal and ferromag- 
net, respectively, A N is the cross section of the normal 
metal and A J is the contact area. For our samples, we 
estimate (K F rs 30) > {K N w 10). In that limit, and 
neglecting contact resistances and inhomogenous current 
distributions, Ai?N is given by 
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(1) 



A joint description of both local and non-local sub-gat 
conductance in the superconducting state is given by |23J 
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where Ia,b, Ua,b and Ga,b are the current, voltage and 
local conductance of the injector (A) and detector (B) 
contact, Gcar and Gec are the non-local conductances 
due to CAR and EC processes, and using Ga, Gb 3> 
Gcar, Gec throughout. Inverting (0 for the case Ib = 
(voltage detection at B) we find the non-local resistance 



R s = U B /h 
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CAR 



G ec )/(GaGb). (3) 



Quantitative predictions for a diffusive superconductor 
yield, to lowest order in transmission coefficients [24| 



with the dirty limit predictions calculated from the re- 
sistivity. Thus, our data in the superconducting state 
are in good quantitative agreement with the predicted 
behavior. 

To conclude, we have investigated experimentally non- 
local sub-gap transport in superconductors-ferromagnet 
hybrid structures, and shown evidence for the superpo- 
sition of crossed Andreev reflection (CAR) and electron 
cotunneling. CAR has been discussed as a sensitive probe 
of superconducting order parameters |2d| . In our case of 
a dirty (hence isotropic) BCS superconductor, the order 
parameter enters only through an exponential decay with 
£. With sufficient spatial resolution an investigation of 
unconventional superconductors appears feasible. More 
recent proposals in the context of quantum information 
processing (see and references therein) have pointed 
out that CAR, if viewed as the decay of a Cooper pair 
into two electrons in different normal metal leads, cre- 
ates a spatially separated entangled electron pair. This 
entanglement can be probed by measuring noise cross- 
correlations in our structures. 

Useful discussions with B. J. van Wees, C. Strunk 
and M. Aprili are gratefully acknowledged. This work 
was supported by the Deutsche Forschungsgemeinschaft 
within the framework of the Center for Functional Nanos- 
tructures. 
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where G A B are the spin-dependent contact conduc- 
tances, and N, D and £ are the density of states at 
the Fermi level, the diffusion coefficient and the dirty 
limit coherence length of aluminum, respectively. Us- 
ing Einstein's relation Ne 2 D = l/p N , we find Rs ~ 
zLPj^Pqttp^ exp( — dj £) / 4d, where + and — are for par- 
allel and antiparallel configuration, and Pa,b = (G A B — 
G ab )/Ga,b are interfacial spin polarizations, and with 
P = Pa = Pb 



Ai? s 



p2 n N exp(~d/£_) 
2 P d 
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Recent predictions for high transparency interfaces in- 
dicate the same distance dependence and an amplitude 
within a factor of three of the low-transparency case • 
Fits of and JSJ to our data in the normal and su- 
perconducting state, respectively, are shown in Fig. 0] 
Within the experimental scatter, the data are described 
well by the fits. We find Ag = 500 - 1000 nm, consistent 
with values for aluminum obtained by spin-valve experi- 
ments [2(|, and a bulk spin polarization p « 20 %. The 
interfacial spin polarization P sa 50 % obtained from the 
fits in the superconducting state is in excellent agree- 
ment with P w 45 % obtained from point contact spec- 
troscopy on iron 0, 0] , and £ « 300 nm is consistent 
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